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Abstract: We study equilibrium droplets in two-phase systems at parameter values corresponding 
to phase coexistence. Specifically, we give a self-contained microscopic derivation of the Gibbs- 
Thomson formula for the deviation of the pressure and the density away from their equilibrium 
values which, according to the interpretation of the classical thermodynamics, appears due to the 
presence of a curved interface. The general — albeit heuristic — reasoning is corroborated by a 
rigorous proof in the case of the two-dimensional Ising lattice gas. 



1. Introduction 

1.1 The problem. 

The description of equilibrium droplets for systems with coexisting phases is one of the out- 
standing achievements of classical thermodynamics. Standard treatments of the subject highlight 
various formulae relating the linear size of the droplet to a specific pressure difference. One of 
these, called the Gibbs-Thomson formula, concerns the difference between the actual pressure 
outside the droplet and the ambient pressure of the system without any droplets. (Or, in the ter- 
minology used in classical textbooks, "above a curved interface" and "above a planar interface," 
respectively.) The standard reasoning behind these formulae is based primarily on macroscopic 
concepts of pressure, surface tension, etc. But, notwithstanding their elegance and simplicity, 
these derivations do not offer much insight into the microscopic aspects of droplet equilibrium. 
The goal of the present paper is to give a self-contained derivation of the Gibbs-Thomson formula 
starting from the first principles of equilibrium statistical mechanics. 

While straightforward on the level of macroscopic thermodynamics, an attempt for a micro- 
scopic theory of droplet equilibrium immediately reveals several technical problems. First of all, 
there is no obvious way — in equilibrium — to discuss finite-sized droplets that are immersed in 
an a priori infinite system. Indeed, the correct setting is the asymptotic behavior of finite sys- 
tems that are scaling to infinity and that contain droplets whose size also scales to infinity (albeit, 
perhaps, at a different rate). Second, a statistical ensemble has to be produced whose typical 
configurations will feature an equilibrium droplet of a given linear size. A natural choice is the 
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canonical ensemble with a tiny fraction of extra particles tuned so that a droplet of a given size 
is induced in the system. A difficulty here concerns the existence of a minimal droplet size as 
will be detailed below. Finally, for the specific problem at hand, the notions of pressure "above 
a curved interface" and "above a planar interface" have to be reformulated in terms of micro- 
scopic quantities which allow for a comparison of the difference between these pressures and the 
droplet size. 

Some of these issues have previously been addressed by the present authors. Specifically, 
in [4,5], we studied the droplet formation/dissolution phenomena in the context of the canonical 
ensemble at parameters corresponding to phase coexistence and the particle density sUghtly ex- 
ceeding the ambient limiting rarefied density. It was found that, if V is the volume of the system 
and 5N is the particle excess, droplets form when the ratio (SNY'^^^^^'^ /V is of the order of 
unity. In particular, there exists a dimensionless parameter A, proportional to the thermodynamic 
hmit of this ratio, and a non-trivial critical value Ac, such that, for A < Ac, all of the excess will 
be absorbed into the (Gaussian) fluctuations of the ambient gas, while if A > Ac, a mesoscopic 
droplet will form. Moreover, the droplet will only subsume a fraction Aa < 1 of the excess 
particles. This fraction gets smaller as A decreases to Ac, yet the minimum fraction does not 
vanish. It is emphasized that these minimum sized droplets are a mesoscopic phenomenon: The 
hnear size of the droplet will be proportional to <c V^^'^ and the droplet thus occupies 

a vanishing fraction of the system. Note that the total volume cannot be taken arbitrary large if 
there is to be a fixed-size droplet at all. 

The droplet formation/dissolution phenomena have been the subject of intensive study in last 
few years. The fact that d/{d + 1) is the correct exponent for the scale on which droplets first 
appear was shown rigorously in [15] (see also [21]); a heuristic derivation may go back at least 
to [3]. The existence of a sharp minimal droplet size on the scale FV(<i+i) 

was described in [22], 

more recently in [4, 25] and yet again in [2]. In the context of the 2D Ising system, a rigorous 
justification of the theory outlined in [4] was provided in [5]. We note that the existence of a 
minimal droplet size seems to be ultimately related to the pressure difference "due" to the pres- 
ence of a droplet as expressed by the Gibbs-Thomson formula. Indeed, from another perspective 
(which is more or less that of [22, 25]), the formation/dissolution phenomena can be understood 
on the basis of arguments in which the Gibbs-Thomson formula serves as a foundation. Finally, 
we remark that although the generation of droplets is an inherently dynamical phenomenon (be- 
yond the reach of current methods) it is possible that, on hmited temporal and spatial scales, the 
equilibrium asymptotics is of direct relevance. 

The remainder of this paper is organized as follows. In the next subsection (Section 1.2) we 
will present an autonomous derivation of the Gibbs-Thomson formula based on first principles 
of statistical mechanics. Aside from our own (modest) appreciation of this approach. Section 1.2 
is worthwhile in the present context because the rigorous analysis develops precisely along these 
lines. In Section 2, we will restrict our attention to the 2D Ising lattice gas, define explicitly the 
relevant quantities and present our rigorous claims in the form of mathematical theorems. The 
proofs will come in Section 3. 
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1.2 Heuristic derivation. 

Let us consider a two-phase system at parameter values corresponding to phase coexistence. We 
will assume that the two phases are distinguished by their densities and, although the forthcoming 
derivation is completely general, we will refer to the dense phase as liquid and to the rarefied 
phase as gas. Confining the system to a (d > 2)-dimensional volume V, we will consider a 
canonical ensemble at inverse temperature /3 and the number of particles fixed to the value 

N = p^V + ip,-pg)SV. (1.1) 

Here, and pg are the bulk densities of the liquid and gas, respectively, and the particle excess 
is 5N = {p^ — Pg)5V with 5V <^ V. Let wi denote the dimensionless interfacial free energy 
(expressed in multiples of /?~^), which represents the cost of an optimally-shaped droplet of unit 
volume, and let x denote the response function, x = y{{N — (N))'^), which is essentially the 
isothermal compressibility. Then, as has been argued in [4], if the parameter 

is less than a critical value Ac = all of the particle excess will be absorbed by the 

background fluctuations, while, for A > Ac, a fraction of the excess particles will condense 
into a droplet. Moreover, the volume of this droplet will be (in the leading order) XaSV, where 
Aa £ [0, 1] is the maximal solution to the equation 



-X-^/'^ = 2A(1 - A). (1.3) 



d 

Note that Aa^ = 2/(d + 1) as advertised; that is to say, the droplet does not appear gradually. 
Furthermore, as is of interest in certain anisotropic situations where the droplet plays a role of an 
equilibrium crystal, the droplet has a particular shape, known as the Wulff shape, which optimizes 
the overall interfacial free energy for a given volume. 

1.2.1 Gibbs-Thomson I: The density. On the basis of the aforementioned claims, we can already 
state a version of the Gibbs-Thomson formula for the difference of densities "due to the presence 
of a curved interface." Indeed, since the droplet only accounts for a fraction, Aa, of the excess 
particles, the remainder (1 — \a){P£ — Pg)^, of these particles reside in the bulk. Supposing 
that the droplet subsumes only a neghgible fraction of the entire volume, i.e., dV ^ V, the gas 
surrounding the droplet wiU thus have the density 

Pg = Pg + (1 - Aa)(p, - Pg)^ (1 + o(l)) . (1.4) 

Here o(l) is a quantity tending to zero as V tends to infinity while keeping A finite (and A > Ac). 
Invoking (1.2) and (1.3), this is easily converted into 

d—1 xwi 1 , ,^ ,N 
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Thus, the density of the gas surrounding the droplet will exceed the density of the ambient gas 
by a factor inversely-proportional to the Unear size of the droplet. This is (quahtatively) what is 
stated by the Gibbs-Thomson formula. 

In order to make correspondence with physics literature, let us assume that the droplet is 
spherical — which is the case for an isotropic surface tension. Then we have 

wi=paSd(^y~ and XaSV = ^r'^ (1.6) 

where a is the surface tension, Sd is the surface area of a unit sphere in and r is the radius of 
the droplet. Substituting these relations into (1.5), we will get 

p-g = Pg + (d-l)-^i(l + o(l)). (1.7) 

Of course, all three formulas (1.4), (1.5) and (1.7) represent the leading order asymptotic in 1/r. 
Higher-order corrections go beyond the validity of the presented argument. 

Remark 1. We note that equation (1.7) differs from the usual corresponding version of the Gibbs- 
Thomson formula in which the x appearing above is replaced by pg. This is due to the approxi- 
mation J< ~ Pg which is justified only in the ideal-gas limit of the rarefied phase. 

1.2.2 Pressures above curved/planar interfaces. Next we turn our attention to the Gibbs- 
Thomson formula for the pressure. Here we immediately run into a complication; while the 
density is a well-defined object in finite volume, the pressure, by its nature, is a macroscopic 
commodity. Thus, strictly speaking, the pressure should be discussed in the context of thermo- 
dynamic Umits. 

In the present context we need to define the "pressure of the gas surrounding a droplet." In 

order to do so, we will consider two canonical ensembles with the same number of particles 
given by (1.1), in volumes V and V + hV , where ^ V . From the perspective of equilibrium 
thermodynamics, these two situations describe the initial and terminal states of the gas undergoing 
isothermal expansion. Standard statistical-mechanical formulas tell us that the change of the 
relevant thermodynamic potential (the Helmholtz free energy) during this expansion is given as 
the pressure times the difference of the volumes aV. Using Zc{N, V) to denote the canonical 
partition function of A'^ particles in volume V, we thus define the relevant pressure pv by 

11 Zc{pgV + {pe-p,)5V,V + AV) 

" °' zc{p,v + ip,-p,)sv,v) ■ ^'-^^ 

For finite V, aV, etc., the quantity pv still depends on aV. As it turns out, this dependence 
(which we will refrain from making notationally expUcit) will annul in any limit V, aV — ^ oo 
with aV/ dV 0, where dV denotes the boundary of V. However, we must consider a limiting 
procedure for which aV also does not "disturb" the droplet. This is a shghtly dehcate subject 
matter to which we will return shortly. 

Our next goal is to give a mathematical interpretation of the pressure "above a planar interface." 
As it turns out (and as is the standard in all derivations), here the correct choice is to take simply 
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the pressure of the ambient gas phase. (See Remark 3 for further discussions.) Using Zq{i^, V) 
to denote the grand canonical partition function, with n denoting the chemical potential, this 
quantity is defined by the (thermodynamic) Umit 

p^ = ^ lim ^\ogZo{lJ.uV), (1.9) 

fj K— >oo V 

Here we have prepositioned the chemical potential to the transitional value, i.e. /i = /it. By well- 
known arguments, this limit is independent of how V tends to infinity provided dV/V tends to 
zero as F — > oo. 

Since we are ultimately looking for an expression for the difference pv — Poo, instead of (1.9) 
we would rather have an expression that takes a form similar to (1.8). We might try to use the fact 
that log Zq^jLi, V) = Pp^V + 0{dV), but then the boundary term will be much larger than the 
actual Gibbs-Thomson correction. We thus have to develop a more precise representation of the 
grand canonical partition function. For simplicity, we will restrict ourselves to the cases when V 
is a rectangular box, in which case we expect to have 

logZoifiu V) = Pp^V + TwaiiSF + O(F^). (1.10) 

Here Twaii denotes a wall surface tension which depends on the boundary conditions. The error 
term represents the contribution from lower-dimensional facets of V, e.g., edges and comers of V 
in d = 3. Such a representation can be justified using low-temperature expansions, see [8], and/or 
by invoking rapid decay of correlations. Of course, this will be discussed in excruciating detail 
in Section 3 of the present paper. 

Using the representation (1.10), we can now write 

which supposes that both V and V + aF are rectangular volumes. 

Our goal is to Umit aF to the values for which the error term is negUgible compared with 

the anticipated Gibbs-Thomson correction. First, supposing that aF ^ F, we find that the 
difference d{V + aF) - dV is of the order aV/V^I'^. Second, assuming that A from (1.2) is 
finite and exceeding Ac (which is necessary to have any droplet at all), we have 5V ~ V^l . 
These two observations show that the contribution of d{y + hV) — dV to the error term in (1.11) 
is indeed negligible compared with {5V)~^^'^. A similar calculation shows that the the second 
part of the error term, f('^~^^/'^/aF, on the right-hand side of (1.11) is negligible compared 
with {5Vy^l'^ provided that 

aF>F^+^. (1.12) 

It is easy to check — see formula (1.23) — that (1.12) can be satisfied while maintaining aF <C 5V. 
This observation will be essential in the forthcoming developments. 

The formulas (1.8-1.11) can be conveniently subtracted in terms of the probability Py(A?^) 
that, in the grand canonical ensemble, there are exactly N particles in volume F. Exphcitly, 
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denoting 



we get 



Here we have applied (1.12) to simplify the error term. 

1.2.3 Gibbs-Thomson II: The pressure. Now we are in a position to derive the desired Gibbs- 
Thomson formula for the pressure. A principal tool for estimating the ratio of the probabilities 
in (1.14) will be another result of [4] which tells us that, in the limit V ^ oo, 

-logFv{pgV + {p,-pg)6V)=wi{5V)'-^{^l + o{l)), (1.15) 

where $^ is the absolute minimum of the function 

$^(A) = + A(l - A)2 (1.16) 

on [0, 1]. Since pgV + {p^ - pg)dV = pg{V + aV) + {p^ - pg)(a SV), where 

a = l-— (1.17) 

we also have, again in the limit V ^ oo, 

-logFv+^v{PsV+iPt - Ps)SV) = wi{aSV)'^{^l^^^+oil)), (1.18) 

where we have introduced the shorthand A (a) = q;~s~ A. 
Supposing that aV ^ 5V, we can write 

^l(a) = - - Aa)' + 0{AV/5V) (1.19) 

Pi — Pg ov 



and thus, to the leading order in aV/ dV, 

PijPv - Poo) = Wl- 



d 1$^ + 1±1a(1-Aa)^ + o(1) 



pg {svy'" [ d "-^^ d '^"^ ' ^^^^ • ^^-^^^ 

After some manipulations involving (1.16) and (1.3), the square bracket on the right-hand side 
turns out to equal ^^^^^^ + o(l). Thus we finally derive 

P{Pv^P^) = —, ^— — -y^l+ol . (1.21) 

d Pe - Pg {XaSV)^^ ^ 

In the case of an isotropic surface tension, formula (1.21) again reduces to 

Pv-Poo = {d- l)^^i (1 + o(l)) . (1.22) 

This is the (leading order) Gibbs-Thomson correction; the one which is usually derived [24, 29] 
by invoking thermodynamic considerations. We note that here the gas-density pg in the numerator 
is fully justified, cf Remark 1 . 
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Remark 2. We note that higher orders in 1/ r — as predicted by the "exponential" Gibbs-Thomson 
formula in classical thermodynamics — go beyond the validity of the formulas (1.15) and (1.18). 
In fact, as a closer look at the F-dependence of 5V and dV suggests, these corrections may 
depend on the choice of the volumes V and V + and on the boundary condition. We further 
remark that both formulas (1.5) and (1.21) have been derived for the situation when a droplet of 
the dense phase forms inside the low-density phase. However, a completely analogous derivation 
works for a droplet of a low-density phase immersed in a high-density environment (e.g., vapor 
bubbles in water). 

Remark 3. Once we have derived the Gibbs-Thomson formula (1.21), we can also justify our 
choice of for the pressure "above a planar interface." First let us note that, in (1.21), can 
be viewed as a convenient normalization constant — subtracting (1.21) for two different volumes, 
say Vi and V2, the quantity completely factors out. Moreover, if Vi <C V2, the contribution of 
the droplet in V2 to such a difference will be neghgible. Thus, in the hmit when V2 ^ 00 and Vi 
stays fixed, pv^ —PV2 tends to py-^^ —p^ as expressed in (1.21). Since also the droplet in V2 becomes 
more and more flat in this limit, p^ indeed represents the pressure "above a planar interface." 

This concludes our heuristic derivation of the Gibbs-Thomson formula. We reiterate that all of 
the above only makes good sense when has been chosen such that 

yi-i+dii <AF<(5y~ W'^TT. (1.23) 

As is easily checked, these inequalities represent a non-trivial interval of values of aV. In the 
next sections, where we will rigorously treat the case of the two-dimensional Ising lattice gas, the 
inequality on the right-hand side will be guaranteed by taking aF = -qSV and then performing 
the limits V 00 followed by 0. 

2. Rigorous results 

2.1 The model. 

Throughout the remainder of this paper, we will focus our attention on the two-dimensional Ising 
lattice gas. The latter refers to a system where each site of the square lattice I? can be either va- 
cant or occupied by one particle. The state of each site is characterized by means of an occupation 
number nx which is zero for a vacant site and one for an occupied site. The formal Hamiltonian 
of the system can be written as 

= -^n^riy- n^rix- (2.1) 

{x,y) X 

Here (x, y) denotes a nearest-neighbor pair on I? and /x plays the role of a chemical potential. 
Note that the Hamiltonian describes particles with a hard-core repulsion and short-range attraction 
(with coupling constant set to unity). 

The Gibbs measure (or Gibbs state) on particle configurations in a finite volume A c is 
defined using the finite- volume version of (2.1) and a boundary condition on the boundary of A. 
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Explicitly, let dA be the set of sites in \ A that have a bond into A and let be the restriction 
of obtained by considering only pairs {x, j/} n A 7^ in the first sum in (2.1) and sites x G A 
in the second sum. If ua G {0, 1}^ is a configuration in A and ng\ is a boundary condition (i.e., 
a configuration on the boundary dA of A), and if =^(nA|^^aA) is the Hamiltonian for these two 
configurations, then the probability of nj^ in the corresponding Gibbs measure is given by 

p-/3^(nA|nsA) 

pnaA,/3,M( ) ^ ^ (2.2) 

^ ^ Z^^-'^{^^,A) 

Here, as usual, /? > is the inverse temperature and the normalization constant, Zq^^'^ {fi, A), 
is the grand canonical partition function in A corresponding to the boundary condition uqa- We 
recall that, according to the standard DLR-scheme [16], the system is at phase coexistence if 
(depending on the boundary conditions and/or the sequence of volumes) there is more than one 
infinite-volume limit of the measures in (2.2). Of particular interest will be the measure in L x L 
rectangular volume Al c 1? and vacant (i.e., n^A^ = 0) boundary condition. In this case we 
will denote the object from (2.2) by p°'^''^. 

As is well known, the lattice gas model (2.1) is equivalent to the Ising magnet with the (formal) 
Hamiltonian 

= - J ^ a^ay -h'^ax, (2.3) 

{x,y} ^ 

coupling constant J = 1/4, external field h = jj. — 2 and the Ising spins (a^) related to the 
occupation variables (rix) via ax = 2nx — 1. The ±-symmetry of the Ising model also allows us 
to identify the regions of phase coexistence of the lattice gas model defined by (2.1): There is a 
value /3c = 2 log(l + \/2) of the inverse temperature such that for (3 > Pc and /x = /xt = 2, there 
exist two distinct translation-invariant, extremal, ergodic, infinite-volume Gibbs states for the 
Hamiltonian (2.1) — a "liquid" state characterized by an abundance of particles over vacancies 
and a "gaseous" state, characterized by an abundance of vacancies over occupied sites. In the 
Ising-spin language, these states correspond to the plus and minus states which in the lattice gas 
language translate to the states generated by the fully occupied or vacant boundary conditions. We 
will use (— )^ and (— )* to denote the expectation with respect to the (infinite-volume) "gaseous" 
and "Uquid" state, respectively. 

In order to discuss the Gibbs-Thomson formula in this model, we need to introduce the relevant 
quantities. Assuming /x = /Xt and (3 > (3^, we will begin by defining the gas and liquid densities: 

P$ = P${P) = {no)p and p^ = p^{p) = {no)p, (2.4) 

where no refers to the occupation variable at the origin. Note that, by the plus-minus Ising 
symmetry, (no)^ = (1 — "-o)* and thus p( + pg = 1. Next we will introduce the quantity x which 
is related to isothermal compressibility: 

X = ^ {{nonx)p - pI) . (2.5) 



A MICROSCOPIC THEORY OF GIBBS-THOMSON FORMULA 



9 



The sum converges for all /9 > /3c by the exponential decay of truncated particle-particle corre- 
lations, \{nxny)°^ — Pg\ < e~l^~^l/^, where ^ = ^(/3) < oo denotes the correlation length. The 
latter was proved in [12,28] in the context of the 2D Ising model. 

The last object we need to bring into play is the surface tension or the interfacial free energy. 
In the 2D Ising model, one can use several equivalent definitions. Since we will not need any of 
them expUcitly, it suffices if we just summarize the major concepts as formulated, more or less, 
in [14,26]: First, for each /3 > /3c, there is a continuous function rg: {n G M^: \n\ = 1} 
(0, cxd), called the microscopic surface tension. Roughly speaking, ra(n) is the cost per length of 
an interface with normal vector n that separates a "gaseous" and "hquid" region. This allows to 
introduce the so called Wulff junctional that assigns to each rectifiable curve = {}pt) in 
the value 



Here nt is the normal vector to at the point (^t. 

The quantity W^idD) expresses the macroscopic cost of a droplet D with boundary dD. In- 
deed, as has been established in the course of last few years [14, 19-21,26,27], the probability 
in the measure p"'*^'^' that a droplet of "liquid" phase occurs whose shape is "near" that of the 
set D is given, to leading order, by exp{— >^(5D)}. Thus the "most favorable" droplet shape 
is obtained by minimizing WpidD') over all D with a given volume. Using W to denote the 
minimizing set with a unit volume (which can be explicitly constructed [13, 17, 30]), we define 



By well-known properties of the surface tension, we have w\{(i) > once f3 > /3c- We note that, 
as in the heuristic section — see Remark 1 — the customary factor 1//3 is incorporated into in 
our definition of the surface tension. 

Remark 4. For those more familiar with the magnetic terminology, let us pause to identify the 
various quantities in Ising language: First, if m*(/3) is the spontaneous magnetization, then we 
have pgiP) = i(l-m*(/3/4))andp^(^) = i(l+m*(/3/4)). Similarly, if x(/3) denotes the mag- 
netic susceptibility in the Ising spin system, then x(/3) = x(/?/4)/4. Finally, the quantity wi{f3) 
corresponds exactly to the similar quantity for the spin system at a quarter of the inverse temper- 
ature. 

2.2 Known facts. 

Here we will review some of the rigorous results concerning the 2D Ising lattice gas in a finite 
volume and a fixed number of particles. In the language of statistical mechanics, this corresponds 
to the canonical ensemble. The stated theorems are transcribes of the corresponding results 
from [5]. 

Recall our notation p°'/^'^ for the Gibbs state in L x L rectangular box A^, and vacant boundary 
conditions on c^A^. Let (ul) be a sequence of positive numbers tending to infinity in such a way 
that v^^"^ /\Al\ tends to a finite non-zero limit. In addition, suppose that {vl) is such that Pg|A£,| -|- 




(2.6) 



wi{/3) = Wf){dW). 



(2.7) 
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ip£ — Pg)vL is a number from {0, 1, . . . , |A/,|} for all L. For any configuration (n-r) in A/,, let 
denote the total number of particles in A^, i.e., 

Nl= (2.8) 

Our first theorem concerns the large-deviation asymptotic for the random variable N^. The fol- 
lowing is a rigorous version of the claim (1.15), which, more or less, is Theorem 1.1 from [5]. 

Theorem A Let [3 > (5^ and let the sequence (vl) and the quantities pg = pg{(3), p^ = 
PiiP), X = and w\ = wi{f5) be as defined previously. Suppose that the limit 

A = ^^^^ lim (2.9) 

exists with A G (0, oo). Then 

lim -±= \ogP°/^''\NL = PglAil + (p, - pg)vL) = -wi inf $a(A), (2.10) 

L— »oo WVl 0<A<1 

where $a(A) = \/A + A(1 - A)^. 

We proceed by a description of the typical configurations in the conditional measure 

Pf'^'{-\NL = p,\Al\ + {p, - Pg)vL), (2.11) 

which, we note, actually does not depend on the choice of the chemical potential. Our characteri- 
zation will be based on the notion of Peierls' contours: Given a particle configuration, let us place 
a dual bond in the middle of each direct bond connecting an occupied and a vacant site. These 
dual bonds can be connected into self-avoiding polygons by applying an appropriate "rounding 
rule," as discussed in [14] and illustrated in, e.g.. Fig. 1 of [5]. Given a contour 7, let 1^(7) denote 
the set of sites enclosed by 7. In accord with [5], we also let diam7 denote the diameter of the 
set V{'y) in the £2 metric on Z^. If F is a collection of contours, we say that 7 G F is an external 
contour if it is not surrounded by any other contour from F. 

While "small" contours are just natural fluctuations within a given phase, "large" contours 
should somehow be interpreted as droplets. It turns out that the corresponding scales are clearly 
separated with no intermediate contours present in typical configurations. The following is es- 
sentially the content of Theorem 1.2 and Corollary 1.3 from [5]. 

Theorem B Let (3 > (3^ and let the sequence [vl) and the quantities pg = Pg{P), Pi = 
PliP), X = >c{p), and wi = wi{P) be as defined previously. Suppose that the limit in (2.9) 
exists with A G (0, 00) and let Ac = ^(3/2)^/^. There exists a number K = K(P,A) < 00 
such that, for each e > and L — ^ 00, the following holds with probability tending to one in the 
distribution (2.11): 

(1) If A < Ac, then all contours 7 satisfy diam7 < KlogL. 

(2) If A > Ac, then there exists a unique contour 70 with 

Xavl{1 - e) < |F(7o)| < Aa^;l(1 + e) (2.12) 
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and 

Pe^AVLil - e) < ^ rix < p^\/\vl{1 + e), (2.13) 
xGy(7o) 

where Aa is the largest solution to the equation 

4A\/A(1 - A) = 1 (2.14) 
in [0, 1]. Moreover, all the other external contours 7 7^ 70 satisfy diam j < K log L. 

Remark 5. We note that, in the case A = Ac, there is at most one large external contour satis- 
fying the bounds (2.12-2.13), or no contour beyond KlogL at all. The details of what exactly 
happens when A = Ag have not, at present, been quantified — ^presumably, these will depend on 
the asymptotic of the sequence vl- 

Remark 6. One additional piece of information we could add about the contour 70 is that its 
macroscopic shape asymptotically optimizes the Wulff functional, see (2.6-2.7). While the shape 
of the unique large contour plays no essential role in this paper (it appears implicitly in the value 
wi) we note that statements of this sort were the basis of the (microscopic) Wulff construction, 
initiated in [1, 14] for the case of 2D Ising model and percolation. These 2D results were later 
extended in [15, 19-21,26,27]. The techniques developed in these papers have been instrumental 
for the results of [5], which addresses the regime that is "critical" for droplet formation. Recently, 
extensions going beyond two spatial dimensions have also been accompUshed [6, 10, 11]. We refer 
to [7] and [5] for more information on the subject. 



2.3 Gibbs-Thomson formula(s) for 2D Ising lattice gas. 

Now we are finally in a position to state our rigorous version of the Gibbs-Thomson formula for 
the 2D Ising lattice gas. We will begin with the formula for the difference of the densities, which 
is, more or less, an immediate corollary of Theorem B. 

Theorem 2.1 Let (5 > (3^ and let the sequence (vl) and the quantities pg = Pg(/3), P£ = PiiP), 
K = x(/3), and wi = wi{(3) be as defined previously. Let A G (0, 00) be as in (2.9). Suppose 
that A > Ac = ^(3/2)^/^ and let Aa be the largest solution of the equation (2.14) in the 
interval [0, 1]. Let Ac^l be the set of configurations {nx)x£AL ^^^^ contain a unique large external 
contour 70 — as described in Theorem B — obeying (2.12-2.13), and whose particle density in the 
exterior ofjo, 

Pext(7o)= ,^ E n., (2.15) 

satisfies the bounds 

9 ~r. - |T./ ^11/2 (1 -^)< Pext(70) -Pg<7i ~ |T/^ Ml/2 (^+^)- (^-^^^ 

^ Pe-Ps 1^(70)1^/^ ^ Pe-Ps 1^(70)1^/^ 

Then, for each e > 0, we have 

lim Pf^''\Ae,L\NL = PsIAlI + {Pi - Pg)vL) = 1. (2.17) 
L— »oo 



12 



M. BISKUP, L. CHAYES AND R. KOTECKY 



Remark 7. We note that, up to the e corrections, (2.16) is exactly (1.5) for d = 2. Indeed, by 
Theorem B we know that |F(7o)| = Xavl{^ + o{l)) and the two formulas are identified by 
noting that 5V corresponds to in our setting. Due to the underlying lattice, the Wulff droplet 
is undoubtedly not circular for any j3 > Pc and the better-known form (1.7) of the (density) 
Gibbs-Thomson formula does not apply. 

In order to state our version of the Gibbs-Thomson formula for the pressure, we will first need 
to define the pressure "above a curved interface" — not to mention the planar interface. We will 
closely follow the heuristic definitions (1.8-1.11). Let us consider a sequence (A^) of squares 
in 7? satisfying 

A'lDAl but A'l^Al (2.18) 
for all L. Let Z^''\n, A) denote the canonical partition function in A with N particles, inverse 
temperature /3 and the vacant boundary condition. This quantity is computed by summing the 
Boltzmann factor, 

expl (3 y n^uyl, (2.19) 



{x,y) 

x,yeA 

over all configurations (ux) with YlxeA '^^ ~ Then we let 

(3 \A', \ A^l °S Z°/{p,\A,\ + ip, - p,)v„ A,y ^ • 

As in the heuristic section, the quantity depends on the sequences (A^), {vl), inverse temper- 
ature f3, and also the boundary condition — all of which is notationally suppressed. 

For the pressure "above a planar interface," again we will simply use the pressure of the pure 
(gaseous) phase. If A c is a finite set, we let Zq^(ij,, A) denote the grand canonical partition 
function in A corresponding to the chemical potential n and vacant boundary condition. Recalling 
that Hi = 2, we define 

p^ = ]- lim log Z°'^(/Xt, Al), (2.21) 
P L->-oo I Ail 

where the limit exists by standard subadditivity arguments. 

Suppose that A > Ag and let us consider the event Be,L collecting all configurations in A^, 
that have a unique "large" contour 70, as described in Theorem B, such that, in addition to (2.12- 
2.13), the volume F(7o) satisfies the inequaUties 

I ^^^7^172(1 - ^) < PiP^-P^) < I iT/r\|i/2 (^ + (2.22) 

^ Pe- Pg l^(7o)r^ ^ Pe- Pg iV'tTojr^ 

Somewhat informally, the event represents the configurations for which the Gibbs-Thomson 

formula for pressure holds up to an e error. The next theorem shows that, as L ^ 00, these 

configurations exhaust all of the conditional measure (2.11): 

Theorem 2.2 Let P > (3c and let the sequence {vl) and the quantities pg = Pg{(3), p£ = 
Pj^{(3), X = x{(3), and wi = wi{f3) be as defined previously. Let A G (0, 00) be as in (2.9). 
Suppose that A > Ac = ^(3/2)^/^ and let Aa be the largest solution to (2.14) in [0, 1]. For 
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each e > 0, there exists a number r/o > such that if (A^) is a sequence of squares in 1? 
satisfying (2.18) and 

lim ^^^^i^ = and lim \^^^\M = ^ e (o, ^o], (2.23) 

L— >oo \ Ai,\ L— >oo Vl 

then 

lim P"'^''^' {B,,l\Nl = PsIAlI + {Pe - Ps)vl) = 1- (2.24) 

Remark 8. As before, since |F(7o)| = \avl{^ + o(l)), the equality (2.24) is a rigorous version 
of (1.21) for the case at hand. The rate at which the limit in (2.24) is achieved depends — among 
other things — on the rate of the convergence in (2.23). We note that the constraints (2.23) corre- 
spond to the bounds in (1.23). In particular, there is a non-trivial set of sequences (A^) for which 
both limits in (2.23) are exactly as prescribed. Finally, the restriction that rj > in (2.23) is due 
to the fact that from [5] we have essentially no control on the rate of convergence in (2. 10). Thus, 
to allow the second limit in (2.23) to be zero, we would have to do a little extra work in order to 
clarify the rate at which the limits in (2.23) and (2.10) are achieved. 

3. Proofs of main results 
3,1 Proofs of Theorems 2.1 and 2.2. 

In this section we provide the proofs of our main results. We will commence with Theorem 2.1: 

Proof of Theorem 2.1. The proof closely follows the heuristic calculation from Section 1.2. Fix 
an e > and let us restrict our attention to particle configurations containing a unique external 
contour 70 and satisfying the bounds (2.12-2.13). Recall the definition (2.8) of the quantity Nl. 
We will show that, under the condition 

NL = Pg\AL\ + {pe-Ps)vL, (3.1) 
any such configuration is, for a suitable e' > 0, contained in Ae',L for all L. Let 

A^ext(7o)= Yl (^-^^ 
a;eAi\F(7o) 

The inequalities in (2.13) then directly imply 

|A?ext(7o) - {Nl - Pe>^AVL)\ < ep^X^VL- (3.3) 

Since we work with a measure conditioned on the event (3.1), we can write 

Nl - Pe^AVL = PgdAil - Xavl) + {pe - /9g)(l - Xa)vl- (3.4) 

But |Al| — Xavl = |Al \ F(7o)| + (|^(7o)| — Xavl) and by (2.12), the second term is no larger 
than eXAVL- Combining the previous estimates, we derive the bound 

|iVext(7o) - PglAi \ F(7o)| - {Pi - Pg)(l - Xa)vl\ < eXAVL, (3.5) 
where we also used (inessentially) that + pg = 1 (and thus pg < 1). 
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The first two terms in the absolute value on the left-hand side represent the difference be- 
tween Pext(7o) and pg while the third term is exactly the Gibbs-Thomson correction. Indeed, di- 
viding (3.5) by I AL\y (70) I and noting that, by definition, A^ext(7o) = Pext(7o) I Al\F(7o) |, we get 



Pext(7o) - Pg- [Pe - Pg, 



\Al\V{jo)\ - \Al\V{^o)\' ^ ■ ^ 

Since both the Gibbs-Thomson correction — which arises from the last term in the above absolute 
value — and the error term on the right-hand side are proportional to ul/IA/, \ 1/(70) |, the desired 
bound (2.16) will follow with some e > once we show that 

. ^ (1 - \a)vl 1 xwi 1 ^ . 

\Al\V{jo)\ ^ Pe-pg VXavl 
To prove (3.7), we note that \Al \ l/(7o)|/| A^l = 1-1- o(l), which using (2.9) allows us to write 

„ , =-^^^(1 + 0(1)), L^oo. (3.8) 

|Al\V^(7o)| {pi-pg)^V^^ " 

Using (2.14) in the form A(l — Aa) = 1/(4V^a)5 we get rid of the factor of A, whereby 
(3.7) follows. Since the o(l) term in (3.7) is uniformly small for all configurations satisfying 
(2.12-2.13), the bounds (2.16) hold once L is sufficiently large. □ 

In order to prove our Gibbs-Thomson formula for the pressure, we will need the following 
representation of the grand canonical partition function: 

Theorem 3.1 Let (3 > (3^ and let be as in (2.21). There exists a number G M and, for 
each 9 G (1, 00), also a constant C{P,9) < 00 such that 

\logZ°/{f,„A) - f3pJA\ - CiilMII < C{p,e) (3.9) 

holds for all rectangular volumes A C 1? whose aspect ratio lies in the interval {6~^, 6). 

Clearly, Theorem 3.1 is a rigorous version of the formula (1.10). Such things are well known in 
the context of low-temperature expansions, see, e.g., [8]. Here we are using expansion techniques 
in conjunction with correlation inequalities to get the claim "down to /3c." However, the full 
argument would detract from the main fine of thought, so the proof is postponed to Section 3.2. 

Proof of Theorem 2.2. We will again closely follow the heuristic derivation from Section 1.2. 
First we note that, using Theorem 3.1, we have 



1 ^ ,^g'^(//.,A'J 



I^Poo - I V , log 



\dA'^\-\dAL\ 2C{P,9) 

-l^-^'l K\A,| +KW ^'-^'^ 



Introducing the shorthand 

PA(iV) = P°''^''^'(^n, = iv), (3.11) 
xeA 

invoking the assumption on the left of (2.23) and applying (2.20), this allows us to write 

1 ^A'^^jpgl^ll + jpl - Pg)vL) ^ ^,^-\l2 

|Ai, W °^ FAjPg|AL| + (p,-PgK) 
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Now, by Theorem A we have 

logPAjpg|Ai| + (p^-pgK) =-u'i($X + o(l))V^^, L^cx), (3.13) 

where is the absolute minimum of $a(A) for A G [0, 1]. As to the corresponding probability 
for A^, we first note that 

PbIAlI + (Pe - Pg)vL = Pg|Ai,| + (p^ - pg)aLVL, (3.14) 

where 

a, = l ^mM. (3.15) 

Pe - Pg VL 

By our assumption on the right-hand side of (2.23), converges to a number a given by a = 
1 ^^r?. Again using Theorem A, we can write 

logFA^{pg\AL\ + {pi- Pg)vL) = -wi{^l3/2^ + o{l))V^^, L^oo. (3.16) 
A simple calculation — of the kind leading to (1.20) — now shows that 

^K^/^A-^A = ^^ ^ + 0(.V% viO, (3.17) 

^ Pe~ Pg V Aa 

while (2.23) implies that 

Plugging these equations, along with (3.13) and (3.16), into (3.12), we have 

l,l^,-,^) = lJV!^'(l^'M + 0(,)), (3.19) 
^Pl-Pg VXavI y rj ) 

where o(l) denotes a quantity tending to zero as L ^ oo while 0(r\) is a quantity independent 
of L and tending to zero at least as fast as r/ in the limit r/ | 0. Equation (3.19) shows that, once L 
is sufficiently large, a particle configuration satisfying the bounds (2.12) from Theorem B will 
also satisfy the bounds (2.22). The limit (2.24) is then a simple conclusion of Theorem B. □ 

3.2 Representation of the partition function. 

The goal of this section is to prove Theorem 3.1. As already mentioned, we will employ two basic 
techniques: cluster expansion and correlation inequalities. The basic strategy of the proof is as 
follows. First we pick a large negative number \xq < //t and use cluster expansion to establish a 
corresponding representation for the partition function Zq^{ij,q, A^). Then, as a second step, we 
invoke correlation inequalities to prove a similar representation for the ratio of the partition func- 
tions {fiQ, Al) and Z^^ {fit, Al). Essential for the second step will be the GHS inequality 
and the exponential decay of correlations for all (3 > (3c. Combining these two steps, the desired 
representation will be proved. 

Let Poo {p) denote the pressure corresponding to the chemical potential p, which is defined by 
the limit as in (2.21) where p^ is replaced by p. (Throughout this derivation, we will keep (3 fixed 
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and suppress it notationally whenever possible.) The first step in the above strategy can then be 
formulated as follows: 

Lemma 3.2 Let (3 > /3c and let Pooi^jj) be as defined above. For each 9 G (1, oo) and each 
sufficiently large negative fj,o, there exists a number rf (/xo) G M and a constant C\ (/?, /xo , ^) < oo 
such that 

|logZ°''^(Mo,A) -/3Poo(/io)|A| -ri°(/io)|aA|| < Ci(/3, ^uq, (3.20) 
holds for each rectangular volume t^dH? whose aspect ratio lies in the interval (0^^, 6). 

To implement the second step of the proof, we need to study the ratio of the partition functions 
with chemical potentials /Ltt and /Lto- Let A be a finite rectangular volume in I? and let 
denote the expectation with respect to the measure in (2.2) with vacant boundary condition. 
Let N\ = X^j-gA Ux- For any jjlq < fii we then have 

H|^^=r(iVA)f"'dM (3.21) 



and 



P{pM-pM)= (no)°''^''^d/i. (3.22) 

where (— )°'/^'/^ denotes the infinite-volume limit (which we are assured exists) of the state {—)°^^'^- 
(Note that (3.22) is true with any infinite-volume Gibbs state substituted.) Combining (3.21- 
3.22), we thus get 

To derive the desired representation, we need to show that the integrand is proportional to \dA\, 
up to an error which does not depend on A. This estimate is provided in the following lemma: 

Lemma 3.3 Let (3 > (3^ and 9 G (1, oo). There exists a constant C2{P, 9) < oo and a bounded 
function : (— oo, /Xt] [0, oo) such that 

\{NKrf'^ - |A|(no)°'^''^ - |5A|t2°(/x)| < C2{(3,9), /x G (-oo,AXt], (3.24) 
holds for each rectangular volume A C whose aspect ratio lies in the interval {9^^, 9). 

Lemma 3.2 will be proved in Section 3.3 and Lemma 3.3 in Section 3.4. With the two lemmas 
in the hand, the proof of Theorem 3.1 is easily concluded: 

Proof of Theorem 3.L Let 9 G (1, oo) and let A be a rectangular volume whose aspect ratio lies 
in the interval {9~^, 9). Fix /xq to be so large (and negative) that Lemma 3.2 holds and let Qi{l^o) 
denote the quantity in the absolute value in (3.20). For each /x G [/xq, /xt], let Q2{l^) denote the 
quantity inside the absolute value in (3.24). Let us define 

/•Ml 

T"waii = '^r(/^o) + / T2(/x)d/x. (3.25) 
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A simple calculation combining (3.20), (3.24) witli (3.23) then shows that 

log Z°/{^i, A) - f3pM\^\ - rl^m = Qi(/xo) + / (52(M)d/x. (3.26) 

Using (3.20) and (3.24), we easily establish that the absolute value of the quantity on right-hand 
side is no larger than C{P, 9) = Ci {/3, ijlq,0) + (fit - no)C2 (/?, 9). □ 

3.3 Cluster expansion. 

Here we will rewrite the grand canonical partition function in terms of a polymer model, then we 
will collect a few facts from the theory of cluster expansions and assemble them into the proof of 
Lemma 3.2. The substance of this section is very standard — mostly siphoned from [23] — so the 
uninterested reader may wish to consider skipping the entire section on a first reading. 

We begin by defining the polymer model. Given a configuration n\ in A, let us call two distinct 
sites of I? connected if they are nearest-neighbors and are both occupied in the configuration n\. 
A polymer is then defined as a connected component of occupied sites. Two polymers are called 
compatible if their union is not connected. A collection of polymers is called compatible if each 
distinct pair of polymers within the collection is compatible. Clearly, the compatible collections 
of polymers are in one-to-one correspondence with the particle configurations. Finally, let us 
introduce some notation: We write P 7^ P' if the polymers P and P' are not compatible and say 
that the polymer P is in A if P C A. 

Let P be a polymer containing A^(P) sites and occupying both endpoints of E{P) edges in 7?. 
We define the Boltzmann weight of P by the formula 

Q^A^) = e^''^''^^^''^''^ (3.27) 

As is straightforward to verify, the partition function Zq^{ii, A) can be written as 

Z°'VA) = ^n WP)' (3.28) 

where the sum runs over all compatible collections of polymers in A. 

This reformulation of the partition function in the language of compatible polymer configu- 
rations allows us to bring to bear the machinery of cluster expansion. Following [23], the next 
key step is a definition of a cluster, generically denoted by C, by which we wiU mean a finite 
non-empty collection of polymers that is connected when viewed as a graph with vertices labeled 
by polymers P G C and edges connecting pairs of incompatible polymers. (Thus, if C contains 
but a single polymer it is automatically a cluster. If C contains more than one polymer, then any 
non-trivial division of C into two disjoint subsets has some incompatibihty between some pair 
chosen one from each of the subsets.) In accord with [23], a cluster C is incompatible with a 
polymer P, expressed hy C^/^ P, if C U {P} is a cluster. 
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In order to use this expansion, we need to verify the convergence criterion from [23]. In present 
context this reads as follows: For some k>0 and any polymer P, 

E CAM(P')e('+'^^^('') < iV(P). (3.29) 

P': 

Since C/3,/u(P) < e'^''''"^^)^'^'^) is true, this obviously holds if is sufficiently large and negative. 
The main result of [23] then says that each cluster C can be given a weight Cp,n{C) (which is 
defined less implicitly in [23]), such that for all finite volumes A c we have 

logZ°'^(/x,A)= C/3,m(Q, (3.30) 

where denotes the set of all clusters arising from polymers in A. Moreover, this expansion is 
accompanied by the bound 

E |CAM(C)|e''^(C)<iV(P), (3.31) 

C: C/P 

where N[C) denotes the sum of A^(P') over all P' constituting C. With (3.30-3.31) in hand, we 
are now ready to prove the first part of the representation of A): 

Proof of Lemma 3.2. First, we will introduce a convenient resummation of (3.30). For each 
polymer P, let ^(P) be the set of sites constituting P. Similarly, for each cluster C, let ^(C) 
be the union of ^(P) over all P constituting C. For each finite ^ C Z^, we let 

^0,i_,iA)= Yl C/3,m(C). (3.32) 

C: J^{C)=A 

Clearly, the weights i?^^^ are invariant with respect to lattice translations and rotations, having 
inherited this property from <^^^^. Moreover, as is easily checked, = unless A is a 

connected set. The new weights allow us to rewrite (3.30) and (3.31) in the following form: 



with 



logZ°'^(/i,A)= Yl (3-33) 

A: AcA 



E \^0,M)\<e-^'' (3-34) 

A: OGA 

\A\>n 

for each n > 0. Here denotes the number of sites in A. 

Now we are in a position to identify the relevant quantities. First, the limiting version of the 
expression (3.33) suggests that the pressure should be given by the formula 

/3pM= Y Tm'^pA^)- (3-35) 

A:OeA 

To define the constant {h) representing the wall surface tension, let EI denote the upper half- 
plane in 7?, i.e., H = {(xi, 2:2) G Z^ : X2 > 0}, and let L be the "line" in 7? corresponding to 
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the boundary of H, i.e., L = {{xi,X2) G : X2 = 0}. Then we define 

A:OeA III I 

Clearly, in order to contribute to r°(/x), the set A would have to have both ^ n H and ^ \ H 
nonempty. On the basis of (3.34) it can be shown that the sums in (3.35) and (3.36) converge 
once (3.29) holds with a k > 0. 

Combining (3.33) with (3.35), we can now write that 

x£A A: xeA ' ' 
AcA 



(3.37) 



xeA A: xeA 
A^A 



Using the fact that A is a connected set and thus AnA7^0andA\A7^0 imply that A n 7^ 0, 
the second term on the right-hand side can further be written as 

^ \Ar\k\ ... V- l^nAl ^p,^{A) 

1^1 A. 1^1 \Ar\dh\ 

A:A<tA ' I x^dAA-.x^A III I oo^ 

xedA A: xeA III \ i 

Here H^; denotes the half -plane in that contains A and whose boundary = dM^ includes 
the portion of the boundary dA that contains x. (Remember that A is a rectangular set and thus 
its boundary dA splits into four disjoint subsets — the sides of A.) 

Let Qi {ji) denote the (complicated) second term on the right-hand side of (3.38). Let ^ be the 
collection of all finite connected sets A<z7?. Notice that, whenever a set ^ G ,2/ intersects dK 
in only one of its sides and A n dh. = Ar\'Lx, then also Ar\ K = Ar\ Hx, and the corresponding 
term in (3.38) vanishes. It follows that, in order for the set A to contribute to the x-th term 
of Qi{ijl), it must contain at least as many sites as is the £°°-distance from x to the sides of dK 
not containing x. Thus, for a given x G dk, a set ^ C 1? can only contribute to Q\ {ji) if A e ,s/ 
and \A\ > dist(x, M \ L^;). 

Since |yl n A|, |^ n ]H[a;| < \A\ and \A n dA\, \A n Lj;| > 1 for any A contributing to Qiin), 
we can use (3.34) to get the bound 

\Qi{f^)\ <Y1 E ^ E e-«d-t(-'^^\^^). (3.39) 

xedA AGS!/, xeA xedA 
\A\>dist{x,dA\I.a:) 

Choosing k > 0, letting G{k) = Yl'^=i < 0°' using Li, L2 G [9~'^L, 0L] to denote the 
lengths of the sides of OA, we can bound the right hand side by 8G{k) + 2Lie'~'^^^ + 2L2e~'^^^, 
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yielding < 8G(K) + 4^Le e^. This in turn can be bounded uniformly in L by a constant 

that depends only on and we thus get the claim of Lemma 3.2. □ 

3.4 Correlation bounds. 

This section will be spent on proving Lemma 3.3. We begin by recalling the relevant correlation 
bounds. Let us extend our notation {—)°if'^ for the expectation with respect to the Gibbs measure 
in A also to the cases when A is not necessarily finite. (It turns out that, by FKG monotonicity, 
such a state is uniquely defined as a Umit of finite-volume Gibbs states along any sequence of 
finite volumes increasing to A.) We will use the notation 

{nx, ny)f^^ = {nxuyy/''' - ''^(n,)f ''^ (3.40) 

for the truncated correlation function. This correlation function has the following properties: 

(1) For each n < n' < Ht and A c A', and all x,y e 1?, 

{nx;ny)f'>'<{nx;ny)l?'>''. (3.41) 

(2) For each P > Pc there exists a ^ = < oo such that 

0<K;n,)f'^<e-l^-^l/« (3.42) 

for all jj, < fit, all A c and all x, y G Z^. Here \x — y\ denotes the l^o distance 
between x and y. 

Both (1) and (2) are reformulations of well-known properties of the truncated correlation func- 
tions for Ising spins. Namely, (1) is a simple consequence of the GHS inequality [18], while (2) is 
a consequence of (1) and the fact that the infinite- volume truncated correlation function at /x = /Xt 
decays exponentially once /3 > Pc- The latter was in turn proved in [12,21]. 

A simple consequence of the above observations is the following lemma: 

Lemma 3.4 Let P > Pc- Then there exist constants ai = ai(/3) G (0, oo) and a2 = a2{P) € 
(0, oo) such that 

< {rixYjf''' - {uxy/'" < aie-"2 dist(a;,A'xA) (3 43^ 

holds for all jx < /Xt, all (not necessarily finite) volumes A C A' C and all x E A. 

Proof. See, e.g., formula (2.2.6) from [21]; the original derivation goes back to [9]. □ 
Now we can start proving Lemma 3.3: 

Proof of Lemma 3.3. We begin by a definition of the quantity (/x) . Let H be the upper half-plane 
in 1?, see Section 3.3. Then we define 

r2°(/^) = Y.^{n^0A)w^' - (^o)z^''^), (3.44) 

where (a; 1, 0^2) is a notation for a generic point in 1?. By Lemma 3.4, the sum converges with 
a /^-independent rate (of course, provided n < Ht). 
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Let A be a rectangular volume in 1? with aspect ratio in the interval {9^^, 9). Let us cyclically 
label the sides of A by numbers 1, ... ,4, and define Hi, . . . ,H4 to be the half-planes in 7? 
containing A and sharing the respective part of the boundary with A. Let us partition the sites 
of A into four sets Ai , . . . , A4 according to which Hj the site is closest to. We resolve the cases 
of a tie by choosing the with the lowest j. Now we can write 

{N^)f''-\K\{n,)l^'^ 

= E E ((-^)a''' - (-^)Hf) + E E ((-^)h^ - ■ (3.45) 

j=l xEAj j=l x&\j 

If it were not for the restriction x G A, the second term on the right-hand side would have the 
structure needed to apply (3.44). To fix this problem, let Sj, with j = 1, . . . , 4, denote the half- 
infinite slab obtained as the intersection IHIj_i n Hj n Hj+i, where it is understood that Hq = IHI4 
and H5 = Hi. Clearly, Kj c for all j = 1, ... , 4. Then we have 

E E ((-^)Hf - (-o>^f 

j=i xeAj 

= r^{n)\dA\ - E ii^-y4'' - (^o)°f ''^). (3.46) 

It remains to show that both the first term on the right-hand side of (3.45) and the second term on 
the right-hand side of (3.46) are bounded by a constant independent of n and A with the above 
properties. As to the first term, we note that, by Lemma 3.4, 

\{nx)f'^ - {nxru-'^l ^ «ie-°^<i^^*(-'«AA), (3.47) 

which after summing over x G Aj gives a plain constant. Concerning the second contribution 
to the error, we note that {nx)^^''^ — {1^0)^2'^ is again exponentially small in dist(x, \ Hj). 
As a simple argument shows, this makes the sum over a; G S j \ Aj finite uniformly in A with a 
bounded aspect ratio. This concludes the proof. □ 
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